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(N) ■ ABSTRACT 

Context. The Atacama Large Millimeter/submillimeter Array (ALMA) will have the necessary resolution to observe a planetary gap 
Qh . created by a Jupiter-mass planet in a protoplanetary disk. Because it will observe at submillimeter and millimeter wavelengths, grains 

£jj ■ in the size range 10yum to 1 cm are relevant for the thermal emission. For the standard parameters of a T Tauri disk, most grains of 

this size range are weakly coupled to the gas (leading to vertical settling and radial migration) and the common approximation of 
well-mixed gas and dust does not hold. 

Aims. We provide predictions for ALMA observations of planet gaps that account for the specific spatial distribution of dust that 
^ . results from consistent gas+dust dynamics. 

5^ ■ Methods. In a previous work, we ran full 3D, two-fluid Smoothed Particle Hydrodynamics (SPH) simulations of a planet embedded 

' in a gas+dust T Tauri disk for different planet masses and grain sizes. In this work, the resulting dust distributions are passed to 

the Monte Carlo radiative transfer code MCFOST to construct synthetic images in the ALMA wavebands. We then use the ALMA 
simulator to produce images that include thermal and phase noise for a range of angular resolutions, wavelengths, and integration 
times, as well as for different inclinations, declinations and distances. We also produce images which assume that gas and dust are 
well mixed with a gas-to-dust ratio of 100 to compare with previous ALMA predictions, all made under this hypothesis. 
Results. Our findings clearly demonstrate the importance of correctly incorporating the dust dynamics. We show that the gap carved 
■ by a 1 M; planet orbiting at 40 AU is visible with a much higher contrast than the well-mixed assumption would predict. In the case 

of a 5 Mj planet, we clearly see a deficit in dust emission in the inner disk, and point out the risk of interpreting the resulting image 
as that of a transition disk with an inner hole if observed in unfavorable conditions. Planet signatures are fainter in more distant disks 
, but declination or inclination to the line-of-sight have little effect on ALMA's ability to resolve the gaps. 

Conclusions. ALMA has the potential to see signposts of planets in disks of nearby star-forming regions. We present optimized 
^■■^ ■ observing parameters to detect them in the case of 1 and 5 Mj planets on 40 AU orbits. 

Key words. Protoplanetary disks - Planet-disk interactions - Methods: numerical - Submillimeter: planetary systems 
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1. Introduction scenario, giant planets form in two steps. A rocky core is first 

formed by progressive accumulation of planetesimals until a 

The study of exoplanets has steadily gained momentum in the critical mass ( „ 10 Me ) is reached. Then rapid capture of nebular 

past two decades, in particular since the confirmed identifier gas can occuri - m a mnaway fashiolli until the planet reaches its 

tion of the first exoplanet around the solar-like star 51 Peg final mass In the disk instability scenario, a gravitationally un- 

( |Mayor&Queloz|1995|) . More than 700 planets located in -500 stable disk f ragmen ts into self-gravitating clumps which then di- 

planetary systems are now confirmed^ A direct consequence rectly form gas giant planetfj The relevance of each mec hanism 

of these surveys is that more reliable estimations of the sta- to explain the current observa tions of exoplanets, especially at 

tistical properties of planetary systems are becoming available large orbkal radiii depends largely on th e exact propertie s of the 

(e.g., distributions of orbital distances and planet masses, see underlying prot oplanetary disk (see e.g. lMordasini et~all20T2h . 
Schneide r et al.H201 II) . 

Two scenarios currently co mpete for giant gase ous planet Although circumstellar disks of gas and dust were imaged 

formation: core-accretion ( e.g. lAlibert et al.l 120051) and disk direct l y roughly at the same time as the exoplanets (iDutrev et al.l 

gravitational instability (e.g. lBossll20lU and references therein), Il994t iBurrows et al.l 1 19961) it is only recently that the link 

each with its own merits and limitations. In the core-accretion between them and the exoplanets was confirmed directly by 

observations when planet s and disks were ima ged simultane- 

1 http://exoplanet.eu ously around /? Pictoris dLa grange et alJ l2010l) and HR 8799 



1 



J.-F. Gonzalez et al.: Planet gaps in the dust layer of 3D protoplanetary disks. II 



dMarois et al ] 120081 l2010h . Both systems are however rather 
evolved and gas giant formation has likely ceased in these now 
gas-poor debris disks. Interestingly, the case of HR 8799 with 
several planets orbiting between 14 and 68 AU suggests that the 
formati on of giant p l anets on wide orbits is a strong possibility. 
Indeed. lOuanz et al] (1201 2|) recently found that a larger fraction 
of stars than previously assume d may harbor gas giant planets 
at larger orbital separations and lLambrechts & Johansenl d2012l) 
showed that such planets can be formed by the core accretion 
mechanism well within the disk lifetime. 

ALMA, when completed, will have the capacity to observe 
these wide planetary systems at a younger age and explore a 
parameter space that is complementary to that probed by ra- 
dial velocity and transit techniques. In this paper we show 
that signs of the presence of these planets on wide orbits, and 
specifically the gaps they carve in the disk, will be detectable 
with ALMA. The process of gap formation depends on the 
planet mass, the disk surface density and the size of d ust grains 
dPaardekooper & Melle"mall2004t iFouchet et aDl2007h and mea- 
suring their width and de pth can provide constraints on the 
underlying disk structure dCrida et alj|2006l) . Previous predic- 
tions indicate that ALMA should be able to detect such gaps 
dWolfetal.ll2002t IWolf & D' Angelc1l2()05l) . Most of these pre- 
dictions were made assuming the very best performance ALMA 
will provide, i.e., longest baselines, shortest wavelengths, no or 
low phase noise. The dust responsible for the emission was also 
assumed to be well-mixed with the gas, irrespective of the dust 
particle size. Here we will use the results of numerical simu- 
lations dFouchet et alj|20iol hereafter Paper I) where the dust 
distributions are calculated for each grain size using a two-phase 
SPH approach to produce ALMA synthetic images. We will as- 
sess the different planet signatures due to different dust distribu- 
tions as a function of the grain size. We will explore the ability of 
ALMA to detect these gaps, but also to measure the disk proper- 
ties both inside and outside the planet's orbit for various combi- 
nations of observing time, angular resolution and wavelength in 
different star forming regions in order to estimate the optimum 
survey strategy. 

2. Methods 

In this section, we present the three numerical tools which we 
use to produce synthetic ALMA images of disks with planetary 
gaps. We first introduce our hydrodynamic code, from which we 
obtain the three-dimensional dust and gas distributions for a pro- 
toplanetary disk with an embedded planet. This dust distribu- 
tion is fed into a radiative transfer code to produce raw intensity 
maps. Finally, to make predictions of what ALMA will see for 
a given wavelength, angular resolution, and integration time for 
disks at a specific distance, declination and inclination, we use 
the ALMA simulator from the CASAEI (Common Astronomy 
Software Applications) package. 

2. 1 . Hydrodynamic simulations 

In a previous work dPaper J) , we ran simulations of a disk of 
mass 0.02 M Q extending from 4 to 160 AU around a solar-mass 
star, typical of Classical T Tauri Stars. These were run with our 
two-phase 3D SPH code (iBarriere-Fouchet et al.ll2005h compris- 
ing gas and dust, in which the dust is treated as a pressureless 
fluid. The initial dust-to-gas ratio is set to 0.01. Both phases in- 
teract via aerodynamic drag in the Epstein regime. Grain sizes 
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Fig. 1. Results of the hydrodynamic simulations: volume density 
maps in midplane (top) and meridian plane (bottom) cuts of the 
disks. The leftmost column shows the gas density and the three 
rightmost columns show the dust density for lOOjurn, 1 mm and 
1 cm grains, from left to right. The rows show simulations with 
M p = 1 and 5 Mj, from top to bottom. The colorbar applies to all 
plots. 



are fixed (grain growth and fragmentation are neglected) and 
self-gravity is not included in this study. We use a vertically 
isothermal temperature profile, T = Tor -1 , with an initial disk 
aspect ratio of H/r = 0.05. The initial density profile is flat, 
Z(r) = So- 

We embed a planet on a circular orbit at 40 AU in a gas-only 
disk, which we evolve for 8 planetary orbits. We then inject the 
dust phase by overlaying dust particles on the gas and evolve the 
system for a total of 104 planetary orbits, following the influ- 
ence of the planet on the vertical settling and radial migration 
of solids. Simulations were ran for different grain sizes: 100 /urn, 
1 mm, 1 cm (smaller grains being well coupled to the gas) and 
planet masses: 0.1, 0.5, 1 and 5 Mj. For more details, see lPaperU 
In this work, we use the simulation results for planet masses of 
1 and 5 Mj. Fig. Q] shows the results of the hydrodynamic simu- 
lations forboth planet masses and all three grain sizes. 

In iPaper J, we showed that the resulting gap is deeper and 
wider for larger grains, as well as for more massive planets. 
The 1 Mj planet only carves a shallow gap in the gas while it 
opens a deep gap in the dust phase. The more massive 5 Mj 
planet carves a deep gap in both the gas and dust phases and 
also traps the 1 cm grains in corotation inside the gap (visible 
as the horseshoe-shaped feature in Fig.Q]), whereas they pile up 
along the spiral wave in the outer disk, causing an asymmet- 
ric structure. Different radial migration efficiencies for different 
grain sizes cause the varying outer disk radius and inner disk dust 
density (in this paper, we call inner and outer disk the disk re- 
gions interior and exterior to the planet gap, resp.). A n extens ive 
discussion of the resulting dynamics can be found in lPaper j . 



2.2. Coupling of the hydro simulations and radiative transfer 
code 

Synthetic images are computed using MCFOST, a 3D contin- 
uum radiative transfer code based on the Monte Carlo method 
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dPinte et al.ll2006l 120091) . The code propagates photon packets 
whose transport is governed by successive scattering, absorption 
and re-emission events determined by the local dust properties. 
The temperature stru cture is calculated using th e immediate re- 
emission algorithm o flBjork man & Woodl(l200l combined with 
a continuou s deposition of energy to estimate the mean intensity 
(lLucvll999l) . We assume dust grains to be homogeneous spheres 
(Mie theory) and adopt t he po rous dust grain optical properties 
from lMathis & Whiffed d!989l model A), with an average grain 
density of 0.5 g crrT 3 . 

To produce synthetic thermal emission maps, we need a pre- 
scription of the three-dimensional dust distribution at every point 
in the disk and thus we need to transform the SPH particle de- 
scription to a grid description. This transformation is directly 
performed on the MCFOST spatial grid. The dust mass in any 
given cell of the radiative transfer code is obtained by summing 
the mass of the neighbouring SPH particles, multiplied by the 
corresponding SPH kernels (see App.lAlfor a discussion of the 
sampling). We use the distributions obtained from our SPH sim- 
ulations for the 100 //m, 1 mm and 1 cm grains. Grains smaller 
than lO O^m, whose spatial distribution was not simulated in 
iPaper J, are strongly coupled to the gas. We therefore assume 
that all grains of 10 //m and smaller are well mixed with the gas 
and follow its spatial distribution. 

Following this procedure, we obtain the spatial density for 
4 grain sizes (lOyum and below, 100 yum, 1 mm and 1 cm). We 
then compute the spatial densities of grains of all other sizes a 
between a m \ n = 0.03 //m and a max = 1 cm by performing a linear 
interpolation in the (log a, log p) plane. The relative amount of 
each grain size is obtained by normalizing the densities so that 
the grain size distribution integrated over the whole disk follows 
dn(a) oc or 3 - 5 da, and the total dust mass is equal to 0.01 of the 
gas mass. We call this setup the "dynamic" case. This procedure 
has been successful ly applied to stu dy dust vertical settling in 
the disk of GG Tau dPinte et al]l2007l) . 

In order to assess the impact of gas and dust decoupling on 
the simulated ALMA images, we also produce grain size dis- 
tributions under the assumption that grains of all sizes are well 
mixed with the gas. In this case, we simply use the resulting SPH 
gas spatial distribution as a proxy to any grain distribution and 
call this setup the "well-mixed" case. 

The inner radius of 4 AU only has a limited impact on the 
disk temperature structure since the disk is optically thick to 
stellar radiation in the radial direction. In the regions probed by 
ALMA (> 15 AU) that we discuss in this work, the heating is 
due to reprocessing of the stellar radiation absorbed in the sur- 
face layers. The temperature structure at large scales is therefore 
not affected by the exact value of the inner radius (see App.lBl). 
As a consequence, only the central beam of the synthetic ALMA 
map depends on our inner radius of 4 AU. 

2.3. Producing synthetic ALMA images 

Images of the disk as it would be observed by ALMA are com- 
puted from the MCFOST maps using the CASA ALMA simula- 
tor (release 3.2, build 15111). The simulator first computes com- 
plex visibilities at each of the (m,v) points sampled by the con- 
sidered ALMA configuration and then adds thermal and phase 
noise before reconstructing an image using the CLEAN algo- 
rithm. Thermal noise is computed by constructing an atmo- 
spheric model from the site characteristics with a given value 
of Precipitable Water Vapor (PWV) and a ground temperature of 
270 K. The atmospheric transparency a bove Chajnantor has been 
monitored over the period 1973-2003 (lOtarola et al.ll2005l) : the 



Table 1. Sky conditions at Chajnantor from lOtarola et al.l 
(l2005h : percentiles of the distribution of the 225 GHz sky op- 
tical depth T225 with its corresponding value and the equivalent 
amount of precipitable water vapor (PWV). 



Percentile 


T 225 


PWV (mm) 


10% 


0.025 


0.30 


25% 


0.0383 


0.58 


50% 


0.0625 


1.08 


75% 


0.1174 


2.22 



quartiles of the distribution of the 225 GHz sky optical depth 
over this period are r 22 5 = 0.0383, 0.0625 and 0.1174. The 
CASA simulator uses a linear relation between T225 and the 
PWV. The corresponding PWVs for the 3 quartiles are 0.58, 
1.08 and 2.22 mm. The best 10% of sky conditions corresponds 
to T225 - 0.025 and a PWV of 0.30 mm. These data are summa- 
rized in TableQ] In the following, we mostly focus on the median 
sky quality and use 1.08 mm of PWV unless otherwise stated. 

Because the troposphere is refractive, atmospheric turbu- 
lence introduces phase delays which are not uniform between 
antennas. This will primarily affect the most extended configura- 
tions and highest frequencies. Two nearby antennas will roughly 
see the same part of the atmosphere and differential phase de- 
lays will be small, whereas well-separated antennas see different 
parts of the atmosphere with uncorrected deformations, result- 
ing in larger phase delays. The larger opacity of water vapor at 
higher frequencies also results in larger phase delays. We use the 
CASA tool settrop to build a 2-dimensional phase screen accord- 
ing to a fractional Brownian motion, allowing us to create real- 
istic correlations in time and space between antennas. In prac- 
tice, this creates a fluctuation of precipitable water vapor screen 
(we assume a relative amplitude of 15%) which is used to cal- 
culate the corresponding phase delay as a function of frequency. 
The phase screen is blown across the array with a wind speed of 
7ms" 1 . ALMA will use Water Vapor Radiometers (WVRs) to 
correct for the atmospheric phase delay in real time during the 
observation. Their efficiency is not yet fully known, and in the 
following we present simulations with and without phase noise. 
We do not include the correction by the WVRs. The quality of 
the ALMA data will lie between these two extreme cases, de- 
pending on the quality of this correction. 

For this work, we have chosen to present reconstructed im- 
ages from simulated ALMA observations rather than visibility 
plots. This is comm on in the literature (e.g. lCossins et al.ll2.010b 
iRegalv et al.l l2012h and therefore easier to compare our results 
with previous predictions of pla net gap observations dWolf et al.l 
I2002t IWolf & D'Angeloll2005h . Images are also a convenient 
way to reveal the instrument signature when the input disk struc- 
ture is known, especially in the favorable case of nearly pole-on 
disks. Finally, the imaging capabilities offered by a large num- 
ber of antennas has been a driver for the building of ALMA. 
However, observers having obtained ALMA data, and without 
any a priori knowledge of the source, will mostly r ely on visibil- 
ities t o understand its structure (similarly to, e.g.. lHughes et al.l 
120071) . 

2.4. Parameter space 

For each of the synthetic maps computed with MCFOST, we 
produce a grid of ALMA images for a range of observational 
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Fig. 2. Raw synthetic images for the 1 Mj planet. Left: well- 
mixed case, center: dynamic case, right: radial brightness pro- 
files of the disks in the well-mixed case (dashed black line) 
and the dynamic case (solid red line). From top to bottom, 
A = 350 yum, 850 fim, 1 .3 mm and 2.7 mm. 

parameters. We choose a reference disk with an inclination 
i = 18.2° (the first in the set of 10 equally spaced values in co- 
sine used in Sect. l3.6t at a declination of -23° and a distance 
of 140pc, and vary the integration time (from lOmin to 8h) 
and angular resolution (from 0.05" to 0.5") for four different 
wavelengths: 350//m, 850//m, 1.3 mm and 2.7 mm. We then de- 
termine the optimal observing mode for all wavelengths (where 
optimal means best compromise between S IN, reasonable inte- 
gration time and angular resolution which recover the features of 
the disk), and then use this observing mode to compare images 
from the well-mixed and dynamic cases at different wavelengths 
for both planet masses. We follow by varying the disk inclina- 
tion and finally, we investigate the effects of source declination 
and distance to determine the observability of planetary gaps in 
different star-forming regions. Table|2]lists the full set of param- 
eters we explored to produce all the ALMA simulated images 
presented in this paper. 

3. Results and discussion 

3. 1 . Raw synthetic images 

In Figs. [2] and [3] we present the raw synthetic images produced 
by MCFOST for both the well-mixed and dynamic cases for four 
different wavelengths and both planet masses. An obvious dif- 
ference between the well-mixed and dynamic results is that the 
extent of the dust disk varies with wavelength in the dynamic 
case but remains fixed for the well-mixed case. The results of 
the hydrodynamic simulations shown in Fig. Q] demonstrate that 
the radial extent of the dust layer varies with grain size, due to 
the varying migration efficiency with grain size. In the synthetic 
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Fig. 3. Same as Fig.|2]for the 5 Mj planet. 



images, all the grain sizes contribute to each wavelength, but 
the grains with sizes closest to the observing wavelength will 
contribute the most. We therefore naturally see the smaller disk 
radius for the longer wavelength in the dynamic case. 

In the case of the 1 Mj planet, the gap is much more visible 
in the dynamic case than in the well-mixed case (Fig.|2]i. This is 
mainly due to the increase in brightness at the gap edges, which 
results in a stronger contrast in the dynamic than in the well- 
mixed case. The hydrodynamic simulations show that this results 
from the pile-up of dust grains in the pressure maxima at the 
gap edges. The brightness contrast increases with wavelength 
in the dynamic case because at longer wavelengths the larger 
grains, which are more effectively cleared out of the gap, are 
the dominant contributors. In the well-mixed case, grains of all 
sizes follow the gas and the resulting contrast is therefore the 
same for all wavelengths. The same contrast effect is seen for 
the 5 Mj planet between the gap and the outer disk (Fig.[3j. For 
the more massive planet, another obvious difference between the 
well-mixed and dynamic models is that features characteristic 
of one particular grain size have different relative contributions 
as the wavelength varies in the dynamic case. The inner disk, 
which only contains lOOpm grains, naturally appears brighter at 
shorter wavelengths. (The larger grains migrate more efficiently 
and hence have almost all been accreted onto the star, see Fig.Q]) 
Similarly, the particles in corotation, which comprise only 1 cm 
grains, appear brighter at longer wavelengths. 

Although both planets produce a spiral wave in the gas 
phase, which is mirrored in the raw images for the well-mixed 
case, it is weaker and barely visible for the 1 Mj planet. In the 
dust phase, grain s accum ulate in the pressure maxima along the 
spiral wave (see iPaper J) , producing correspondingly stronger 
features in the images for the dynamic case. No significant en- 
hancement of the weak spiral is discernible for the 1 Mj planet, 
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Fig. 4. Simulated observations for the 1 Mj planet in the dynamic case at A — 850 fim. From left to right: integration time of 
lOmin, 30min, lh, 2h, 4h and 8h. From top to bottom: angular resolution of 0.05", 0.10", 0.15", 0.30" and 0.50" (the ALMA 
configuration is adapted to match the angular resolution for each wavelength). The scale on each image is in arcseconds, with the 
beam size represented at its bottom left corner, and the colorbar gives the flux in mJy/beam. (Note that the flux scale changes in 
each row, due to the different beam size.) 



but the 5 Mj planet produces a stronger spiral structure in the 
images for the dynamic case. 

While we can see a range of different features in the dynamic 
vs. well-mixed cases, one must remember that these are raw syn- 
thetic images. In order to predict what a real telescope will see, 
we now turn to results from the ALMA simulator. 

3.2. ALMA synthetic images 

We start by investigating the effect of the integration time and an- 
gular resolution on ALMA's ability to resolve features in disks 
with planet gaps. In Fig.|4] we show as an example the ALMA 
simulated images for the 1 Mj planet at A = 850 yum, for expo- 
sure times ranging from lOmin to 8 h and for angular resolutions 
from 0.05" to 0.50". Increasing the integration time increases 
the sensitivity and therefore the signal-to-noise ratio. In partic- 
ular, the background becomes less noisy as the integration time 
increases. However, the gap as well as other features such as 
brightness variation in the disk are already well detectable af- 
ter 1 h and integrating longer does not produce much additional 
detail. An integration time of 1 h thus seems to be the best com- 
promise between being able to detect the essential properties of 



the disk and using a reasonable amount of resources, especially 
in the context of an observatory that is well oversubscribed. 

Figure|4]shows that an angular resolution of 0. 15" or smaller 
is necessary to detect the gap at 850 jum. However, higher an- 
gular resolution requires more extended array configurations, 
with correspondingly decreasing beam sizes, therefore sampling 
smaller areas of the source. The extended array would detect less 
flux, resulting in decreased contrast in the image. We find that a 
resolution of 0.10" gives the best disk-to-gap contrast and adopt 
this value for the rest of this paper. We verify that the optimal 
combination of an integration time of 1 h and angular resolution 
of 0.10" is adequate for all wavelengths studied here, except for 
2.7 mm — see Figs lC. 1 lfC.3l For the longest ALMA waveband, 
the minimum angular resolution required to resolve the gap cor- 
responds to a very extended configuration, and our optimal ob- 
serving parameters result in a very weak signal, making it ex- 
tremely difficult to detect the disk and its features (see App.lCl 
showing figures similar to Fig.|4]for the other wavelengths). 

FigureO shows maps of the signal-to-noise ratio per resolu- 
tion element in the inner and outer disk (i.e. the disk regions 
interior and exterior to the gap, resp.) for both planet masses in 
the wavelength vs. integration time plane for the dynamic case 
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Table 2. Parameters used to produce all ALMA simulated images presented in this paper's figures. 
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Fig. 5. Maps of the signal-to-noise ratio in the inner (top) and outer (bottom) disks in the wavelength-integration time plane for 
the dynamic case, with an angular resolution of 0.10". The signal-to-noise ratio is computed as the median of the flux in a ring 
(between 0.07" and 0.21" for the inner disk and between 0.35" and 0.70" for the outer disk) divided by the rms of the background. 
Left: M p = 1 Mj, right: M p = 5 Mj. Thick dotted, dashed and solid lines show contours for S /N =1,3 and 10. 



with 0.10" angular resolution. These maps, showing very low 
values of S /N for 2.7 mm, confirm our statement of the previous 
paragraph. Only very long integration times would allow the S/N 
to approach or slightly exceed unity at this wavelength, and we 
therefore no longer discuss it in the rest of this paper. For shorter 
wavelengths, Fig. [5] shows that a 1 h integration time is indeed 



enough to reach a S/N larger than 7-8, and even well above 10 
at 850yum (except in the inner disk for M p = 5 Mj, which will 
be discussed in Sect. 13. 4t , ensuring an adequate detection of the 
planetary gaps. 
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Fig. 6. Simulated observation for the 1 Mj planet with an inte- 
gration time of 1 h and an angular resolution of 0.10". Left: 
well-mixed case, center: dynamic case, right: radial brightness 
profiles of the disks in the well-mixed case (dashed black line) 
and the dynamic case (solid red line). From top to bottom, 
A = 350yum, 850jumand 1.3 mm (the 2.7 mm image is not shown 
because it contains only noise.) The scale on the images and ra- 
dial profiles is in arcseconds. 



3.3. "Well-mixed" vs. "dynamic" 

We next compare simulated images with our optimal set of ob- 
serving parameters (1 h integration time, 0.10" angular resolu- 
tion) of the 1 Mj (Fig.[6]) and 5 Mj (Fig.[7j) planet gaps for the 
well-mixed and dynamic cases. Quantitative information is ob- 
tained from the radial brightness profiles of these images, shown 
in the right columns of Figs. [6] and [7] As already seen in the raw 
synthetic images, the disks are more compact in the dynamic 
case, and all the more so for longer wavelengths (down to a ra- 
dius of ~0.8" at 1.3 mm), than in the well-mixed case (with a 
constant radius of ~1.1"). However, a number of features of the 
raw images are lost here. 

In the well-mixed case, the inner disk is too faint to be visible 
for both planet masses at 1.3 mm, showing instead a disk with 
a large inner hole, which is wider for the more massive planet 
(with radii of -0.2" and -OA" for M v = 1 and 5 Mj resp.). At 
850/iin, the 1 Mj planet gap is hard to detect due to the shallow 
(~15%) dip and mild slope of its brightness profile. Its location 
can however be inferred from to the -20% brightness contrast 
between the fainter inner disk and the outer disk. On the other 
hand, the gap is well defined and ~ 1 .9 times fainter than the inner 
disk for the 5 Mj planet. The situation is less favorable at 350 fim 
where the 5 Mj planet gap is deep enough to be detected (~1.5 
times fainter than the inner disk) but its edges are less sharp, 
whereas it is hard to identify a gap in the shallow brightness 
variations for the 1 Mj planet. No additional details can be seen 
to provide constraints on the underlying disk structure. 

In the dynamic case, where the dust features are sharper, the 
planet gaps are very well defined. The 1 Mj planet gap is clearly 
detected at all wavelengths, with sharp transitions to both the in- 
ner and outer disk (only at 350 fim is the transition to the outer 
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Fig. 7. Same as Fig|6]for the 5 Mj planet. 



disk smoother). The disk-to-gap contrast can be estimated from 
the right column of Fig.|6]to be ~ 1 .4 at 350 /an, ~ 1 .6 at 850 fim 
and -2 at 1 .3 mm. For the 5 Mj planet, however, the fainter in- 
ner disk is clearly visible only at 350 fim, whereas it is very faint 
at 850//m (as already seen in the upper right panel of Fig. [5]) 
and barely detectable at 1.3 mm. This, combined with the strong 
contrast between the outer disk and the gap (~2.1, ~2.5 and ~10 
at 350 fim, 850yt/m and 1.3 mm) constitutes a possible source of 
confusion with disks having cleared a large inner hole (see be- 
low). Finally, the corotating 1 cm sized grain population and the 
spiral wave are no longer visible in the ALMA simulated im- 
ages, depriving us of a potential handle on the planet mass or 
angular position. 

These notable differences between images computed in both 
cases constitute one of the main results of our study: the much 
higher contrast of the planet gaps in the more realistic dynamic 
case compared to the well-mixed approximation clearly shows 
that they will be much easier to dete ct than anticipated from 
gas-only si mulations, such as those of Wolf et all d2002l) . This 
was seen bv lPaardekooper & Mellemal (2004 in the case of 2D, 
face-on disks. As we showed in iPaper I , for 3D disks the much 
reduced scale height of the settled dust phase compared to the 
gas further enhances the sharpness of the dust features. As a re- 
sult, they are easier to detect in the simulated images we present 
here. In studying 3D disks it is also possible to investigate the 
effect of disk inclination, see Sect. 13.61 In the rest of this paper, 
we will no longer consider the well-mixed approximation and 
only discuss images obtained from the dynamic case. 



3.4. Gap detectability 

We now assess the detectability of the planet gaps in simulated 
images using our optimal parameters. A single 1-hour ALMA 
image at a well-chosen wavelength (Figs. [6] and [7] show 850//m 
to be a good choice) can be sufficient to detect a planet gap and 
infer the presence of an unseen planet. This is a definite advan- 
tage for future large surveys of star-forming regions. However, 
characterizing the planet which opens the gap will likely re- 
quire multi-wavelength follow-up observations. First, weighing 
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Fig. 8. Simulated observation for the 1 Mj planet with an integra- 
tion time of 1 h and an angular resolution of 0.10" in the dynamic 
case with phase noise included. Left: PWV = 1 .08 mm, center: 
PWV = 0.58mm, right: PWV = 0.30 mm. From top to bottom, 
A = 350/^m, 850 fim, and 1.3 mm. The scale on each image is 
in arcseconds, with the beam size represented at its bottom left 
corner, and the colorbar gives the flux in mJy/beam. 



the planet can not rely solely on the measure of the observed gap 
width, which does not depend strongly on the planet mass (see 
the right columns of Figs.|6]and|7]i, partly because it (as well as 
the gap-to-disk flux ratio) is affected by the beam size. One can 
use in addition the different wavelength dependence of other fea- 
tures produced by planets of different masses, such as the inner 
disk brightness. In our example with 1 and 5 Mj planets, the in- 
ner disk brightness decreases more rapidly towards longer wave- 
lengths for the more massive planet. Furthermore, in the case 
where the inner disk is too faint to be unambiguously detected, it 
is important to be able to differentiate it from a disk with a large 
inner hole. In recent years, a growing number of disks with large 
dus t-free cavities around the central star have been observed (see 
e.g. lAndrews et al 1 1201 lh . The clearing of the inner regions of 
these so-called transition disks is attrib uted either to dynami- 
cal in teractions with a com panion star ( Artvmowicz & Lubowl 
11994 or a massive planet iyarniere et alJl2006h or, more effi- 
ciently, multi-planet systems (iDodson-Robinsqn & Salvkl201 lb. 
or to p hotoevaporative winds (I Alexander et al. 2006; Gorti et all 
l2009b . Our simulated ALMA observations for the 5 Mj planet 



show that the inner disk, while difficult to detect at 850 jum and 
even more so at 1.3 mm, is easily seen at 350/jm (Fig. [7]). Being 
able to image the disk at short wavelengths therefore seems es- 
sential to disentangle disks with planetary gaps from transition 
disks in the dust continuum. Even though the presence of warm 
dust in the inner disk can be hinted at by an IR excess in the 
Spectral Energy Distribution, imaging the inner disk is essential 
to constrain the amount of dust it contains. In the case where 
the inner disk is present but too faint to be unambiguously de- 
tected in reconstructed images, visibilities may show emission 
in the regions close to the star, providing an alternate way of 
discriminating disks with planetary gaps from transition disks. 
A different m e thod f or doing so has recently been proposed by 
ICleeves et al ] (12011b . who claim that the heated inner wall of 



transition disks would produce unique chemical features in the 
gas phase, potentially observable by ALMA. 

It is worth noting that gaps in protoplanetary disks can be 
created without the presence of a planet. Before clearing an 
inner hole, X-ray pho toevaporation ope ns a gap at a few AU 
from solar-type stars dOwen et al 1 12011b but such a feature is 
short-lived (~5% of the disk lifetime, i.e. a few hundred thou- 
sand years) before the whole inner disk is drained on the star. 
In contrast, the planet gaps w e consider here take less than 
20000y r to open (s ee Paper J) and survive until the disk is 
cleared. iPinilla et al.l d2012b considered the effect of a bumpy 
surface density profile and showed that the accumulation of dust 
at local maxima would be seen with ALMA as a succession of 
rings and gaps. They speculate that such a density profile may 
be caused by the magneto-rotational instability, however this re- 
mains to be tested. 



3.5. Effects of phase noise 

The left columns of Figs.[8]and|9]show the effect of phase noise 
on the ALMA simulated images for both planet masses, to be 
compared to the center columns of Figs. [6] and [7] computed with- 
out phase noise. Note that the absolute flux scale is different 
in both sets of figures. As expected, the image deterioration is 
strongest for the shortest wavelengths, to the point where the 
disk is no longer visible at 350//m. Despite some degradation 
at 850yt/m due to the higher noise level, both the disk and the 
gap are still well detected, in fact with a very similar disk-to- 
gap contrast to the images without phase noise. The image at 
1 .3 mm is virtually unaffected by phase noise, and this wave- 
length appears to be the best to recover all the details of the disk. 
Figure[10]shows signal-to-noise maps in the outer disk in the im- 
ages with phase noise for the 1 and 5 Mj planets, to be compared 
with the bottom panels of Fig. [5] Phase noise naturally reduces 
the S IN more strongly as wavelength decreases, whatever the in- 
tegration time. The S IN is almost unchanged at 2.7 mm, where 
it was already very low, slightly reduced at 1.3 mm with values 
~7-9 with an 1 h integration time, degraded at 850 fim to ~3-4, 
and dramatically reduced to values well below unity at 350 /im. 
Note that in the presence of phase noise, integrating longer does 
not improve the image, especially at short wavelengths. In partic- 
ular, Fig.fTOlshows that the S IN at 350/mi even decreases above 
2 h of integration. 

While the detection of planet gaps is still possible at 850 fim 
and good at 1.3 mm, the differential signal-to-noise degradation 
hinders the possibility of using multi-wavelength observations 
to characterize the gap-opening planet. In particular, using the 
lowest wavelength to help separate disks with planet gaps from 
transition disks is no longer possible with the effect of phase 
noise. However, this can be mitigated by the real-time correction 
that will be available thanks to the use of WVRs (see Sect. l2.3l >. 

In practice, it is expected that observations at short wave- 
lengths, and especially at 350/mi, will only be performed in 
very dry weather. We now consider the best 25% and 10% of 
sky conditions and present the corresponding images, computed 
with respectively PWV = 0.58 and 0.30 mm, in the center and 
right columns of Figs.|8]and|9] As mentioned earlier, the effect 
of phase noise is moderate-to-negligible at 850 //m and 1.3 mm, 
respectively, and drier conditions mostly increase the S /N thanks 
to the associated decrease of the thermal noise. At 350 /im, on the 
other hand, even though decreasing the PWV to 0.58 mm barely 
brings the disk out the noise, the improvement is dramatic in the 
driest case considered here, with PWV = 0.30 mm, where both 
the disk and the gap are clearly detected with a contrast similar 
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to the 850 jum image. In the best 10% of sky conditions, both the 
thermal and phase noises are so low that the shortest wavelengths 
are accessible without any necessity for correction and can pro- 
vide the data needed to discriminate transition disks from disks 
with planet gaps even if the WVRs were not efficient. In such 
weather, the use of fully-functional WVRs and fast-switching 
will still improve the S/N. 

3.6. Varying disk inclination 

Since disks are not always viewed pole-on, we now investigate 
the effect of different disk inclinations. The images shown in 
Fig.[TT]are calculated at 10 inclinations equally spaced in cosine, 
i.e. for cos(0 = 0.05, 0.15, ... , 0.85 and 0.95. This corresponds 
to disks randomly oriented in three dimensions, with inclinations 
to the line-of-sight i = 18.2, 31.8, 41.4, 49.5, 56.6, 63.3, 69.5, 
75.5, 81.4 and 87.1°. The planet gap is very easily seen for all 
inclinations up to ; = 75.5°, only hinted at for i = 81.4°, and 
not visible for i = 87.1°. This is because the dust component 
of the disk, responsible for the emission in the continuum, has 
efficiently settled to the midplane and is almost flat (see Fig.[T|). 
There is thus no flared regions in the outer disk that would mask 
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Fig. 11. Simulated observations of disks of varying inclination 
for the 1 Mj planet in the dynamic case at A = 850 //m with an 
integration time of 1 h. The scale on each image is in arcseconds, 
with the beam size represented at its bottom left corner, the disk 
inclination in the upper right corner, and the colorbar gives the 
flux in mJy/beam. 



the gap at high inclinations. Only when the disk is seen nearly 
edge-on does the outer disk hide the inner regions. Gap detection 
is therefore robust with respect to disk inclination. 
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Fig. 12. Simulated observations of disks located in the star- 
forming regions listed in Table [3] for the 1 Mj planet in the dy- 
namic case at A = 850 ^m with an integration time of 1 h. The 
angular resolution for each image is given on the left (it corre- 
sponds to an identical spatial resolution for all star forming re- 
gions), together with the region's name. The scale on each image 
is in arcseconds, with the beam size represented at its bottom left 
corner, and the colorbar gives the flux in mJy/beam. 



3.7. Varying source distance and declination 

ALMA simulated images are almost always shown in the most 
favorable case of a source passing through the zenith at ALMA's 
latitude (as is the case for our reference disk) and for a single dis- 
tance. We present in Fig.[T2]simulated images where we vary the 
declination and distance of the source, taking values for specific 
star-forming regions listed in Table[3] The images are shown for 
the 1 Ms planet at 850 /vm and for a 1 h observing time, without 
phase noise. In order to compare the gap detectability with our 
reference disk at 140 pc observed with an angular resolution of 
0.10", we chose to keep the same corresponding linear resolu- 
tion of 14 AU at the distance to each star-forming region and set 
the angular resolution accordingly. 



star-forming region 


SO 


d (pc) 


Reference 


-23 


140 


Taurus 


+25 


140 


Serpens 


+01 


260 


Ophiuchus 


-24 


120 


Lupus (I) 


-34 


150 


Chamaeleon (I) 


-77 


160 



For the most distant disks in Serpens at 260 pc, the small 
angular resolution required (0.05") is obtained by an extended 
array configuration, resulting in a low signal-to-noise ratio (see 
Sect. l3.21 >. The gap is nonetheless detected, but characterizing 
the planet responsible for it will likely be impractical with these 
observing parameters. The other four star-forming regions are 
approximately twice as close, with similar distances of 120 to 
160 pc, corresponding to required angular resolutions of 0.09- 
0.12". These are close to our optimal value, and naturally pro- 
vide images where the gap is clearly seen. 

Surprisingly, the source's declination has very little influence 
on the gap detectability, even for disks culminating at a rather 
low elevation, like in Taurus (42°) or Chamaeleon (36°). The 
beam is the most elongated for the latter (and as a result samples 
a larger area of the source, making the disk appear brighter) but 
not enough to smear out the sharp gap edges. This very good 
performance compared to disks reaching high elevations is due 
to the excellent uv coverage that ALMA provides even in an ob- 
serving time as short as 1 h and ensures that the disk-to-gap con- 
trast stays high for a large range of source declinations. 

The prospect of detecting planet gaps with ALMA in nearby 
star-forming regions therefore seems excellent. Ophiuchus is ob- 
viously most favorable as it is the closest star-forming region and 
it culminates very close to the zenith at the location of ALMA. 

3.8. Pushing ALMA further 

The detection of planet signatures is a key science driver for 
ALMA, and one can reasonably expect that it will be possible 
to commit a substantial amount of time on promising sources, 
which will be necessary as one wants to reach high angular res- 
olution and high S IN. In the context of trying to obtain the very 
best of what ALMA can offer, we computed simulated ALMA 
images of our reference disk for both planet masses with a 0.05" 
angular resolution, an integration time of 8 h, and no phase noise 
(assuming the WVRs will be able to correct most of it), shown 
in Fig. [13] As expected, the images for both planet masses show 
much more sharply defined gaps, allowing a more accurate mea- 
sure of their width, than for the observing modes discussed 
above. For the 1 Mj planet, the improvement in the disk-to-gap 
contrast is particularly striking at 350 /urn, increasing from ~ 1 .4 
to ~1.8, but no additional features are recovered. On the other 
hand, new details now appear for the 5 Mj planet. At 850 fim 
and 1.3 mm, emission from the corotating grains can be seen but 
unfortunately remains too faint to be used as a constraint on the 
planet mass. However, the brightest parts of the outer disk now 
reveal the spiral wave, especially at 850 jum, which can help con- 
strain the planet's angular position. 

Spiral features have already been detected in visible 
or near-infrared scattered-light images of disks a few hun- 
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Fig. 13. Simulated observation for our reference disk with an in- 
tegration time of 8h and an angular resolution of 0.05". Left: 
M p = 1 Mj, right: M p - 5 M 3 . From top to bottom, A = 350 fan, 
850/im and 1.3 mm. The scale on each image is in arcseconds, 
with the beam size represented at its bottom left corner, and the 
colorbar gives the flux in mJy/beam. 



dred AU from their star, e.g. AB Aur dGradv et alJ Il999t 
iFukagawa etld] |2004|). HP 1 00546 dGradv etal J 120011) and 
HP 141569 A dClampin et al.1 120031) . Subse quent mm obser- 
vations of the well-studied AB Aur disk (Pietuetal 1 120051: 
iLin et al.l l2006l) estimated the Toomre Q parameter well above 
unity, ruling out the gravitational instability scenario that was 
proposed previously and favoring excitation b y a giant planet at 
severa l tens of AU to ~100 AU. More recently. lHashimoto et alJ 
(1201 ll) obtained high-resolution //-band polarized intensity im- 
ages of the AB Aur inner disk and detected two rings separated 
by a gap at 80 AU, possibly caused by a Jupiter-mass planet, and 
concluded that the disk is probably in an early and active phase 
of planet formation. The first exampl e of a direct constr aint from 
observed spiral features comes from lMuto et al.l d2012l) who de- 
tected two spirals in //-band polarized intensity images of the 
HD 135344B disk: using spiral density wave theory, they infer 
the presence of two unseen planets of ~0.5 Mj. ALMA, observ- 
ing at longer wavelengths, will probe structures closer to the 
midplane and bring additional constraints to such studies. 

To push ALMA performances even further, we also com- 
puted simulated images assuming drier conditions giving the 10 
and 25% best values of the sky opacity, corresponding to a PWV 
of 0.30 and 0.58 mm, respectively (see Table[T|). These are shown 
in Fig. [14] for A = 350yt/m only (images at longer wavelengths 
are virtually unaffected). For these lower opacities, the thermal 
noise (as well as the phase noise, not included here) is strongly 
reduced, resulting in images with a much improved S /N at the 
same angular resolution. Figure[l5] shows signal-to-noise maps 
in the inner disk in the images with an angular resolution of 
0.05" and a PWV of 0.30 mm for the 1 and 5 Mj planets. S/N 
values are lower than those in the top panels of Fig. [5] in spite 
of better sky conditions because the smaller beam collects less 
signal in the same amount of time. S/N values above 3 are nev- 
ertheless easily reached in one hour, and Fig. [15] can be used to 
determine the integration time needed to reach a given S /N at a 
given wavelength. 
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Fig. 14. Same as Fig. [13] for A = 350/mi only. Left: M p = 
1 Mj, right: M p = 5 M s . Top: PWV = 0.30mm, bottom: 
PWV = 0.58 mm. 

Pushing ALMA towards its limits therefore seems worth- 
while to better characterize sources of interest. Here, 850//m 
appears to be the most favorable wavelength in average sky 
conditions, recovering the most details in the disk for both 
planet masses. In the driest weather, observations at the short- 
est wavelength of 350 jum are optimal and allow the detection of 
wavelength-dependent features such as the brightness of the in- 
ner disk. However, Fig.[T5lshows that 850/^m remains the wave- 
length giving the best S /N for our disks, even in the best 10% of 
sky conditions. 

4. Conclusion 



In a previous study (see lPaperll) . we have run 3D hydrodynamic 
simulations of the gaps carved in a 0.02 M disk of gas and dust 
by 1 and 5 M s planets at 40 AU from the 1 M central star. We 
followed consistently the dynamics of dust grains in the range of 
sizes contributing the most to the ALMA wavelengths and de- 
scribed the distinctive, sharper features in the dust phase. In this 
paper, we now provide predictions of observations with ALMA 
and examine the detectability of these planetary gaps. We first 
produce raw synthetic images by applying the radiative transfer 
code MCFOST to the hydrodynamic results. We produce two 
kinds of images: the "dynamic" case where we use the true spa- 
tial distribution of dust of different grain sizes obtained from the 
hydrodynamic simulations and the "well-mixed" case where we 
simply use the spatial distribution of gas and assume that the 
dust follows it in order to check the validity of the "well-mixed" 
hypothesis. 

We then produce simulated ALMA images using the CASA 
software package. We investigate the choice of integration time, 
wavelength and angular resolution. We show that an integration 
time of 1 h is enough to ensure a firm detection of our planet gaps 
and that the gain in signal-to-noise as we increase the integration 
time does not warrant the necessary extra observing time on a 
facility that will be heavily oversubscribed. We also show that 
an angular resolution of 0.10" is optimal to resolve the gap with 
the best contrast in most nearby star-forming regions. We then 
show that the use of small or large wavelengths combined with 
these optimal observing parameters give poor to bad results. In 
average sky conditions, we find the wavelength giving the best 
signal-to-noise ratio to be 850 /im when phase noise is neglected 
or 1.3 mm, with slightly lower performance, when it is included. 
However, the use of WVRs should provide real-time correction 
of phase noise and favor 850/mi as the wavelength of choice. In 
the driest weather, phase noise is negligible. 
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Fig. 15. Same as Fig. [5]for the inner disk only, with an angular resolution of 0.05" and a PWV of 0.30 mm. Left: M p = 1 Mj, right: 
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A number of previous stud ies of planet gaps (e.g. lWolf et al.l 

l2002HWolf & D'Angelol2005h . and even recent wor k on Rossby 
vortices in protoplanetary disks dRegalv et al.ll2012T) . have relied 
on hydrodynamic simulations of 2D gas-only disks to produce 
simulations of ALMA observations by reconstructing a three- 
dimensional distribution of dust that follows that of the gas in 
the midplane and is in hydrostatic equilibrium in the vertical di- 
rection before feeding it to a radiative transfer code. In addition 
to neglecting the vertical dust settling, this approach also results 
in smoother dust features, leading to more pessimistic predic- 
tions for the detectability of the disk features. The comparison of 
our simulated images in the "well-mixed" and "dynamic" cases 
show that they are very different and that the more realistic case 
of including a self-consistent treatment of the dust dynamics pro- 
duces much sharper features and more clearly defined gaps. We 
would like to stress here that investigators need to be cautious 
when assessing the observability of a particular structure in their 
source from approximate methods and that a procedure similar 
to our "dynamic" case should be used whenever possible. 

With our "dynamic" approach, we find that the gap carved 
by a 1 Ms planet orbiting 40 AU from its star is easily detected. 
Images for the 5 Mj planet show a high contrast between the 
outer disk and the deeper and wider gap, but the fainter inner disk 
is barely seen with our standard observing parameters. In order 
to detect the regions interior to the gap and remove the possible 
confusion with transition disks with large inner holes, one needs 
to observe at shorter wavelengths in the driest conditions possi- 
ble or use a combination of longer integration time and higher 
angular resolution. In any case, multi-wavelength followup ob- 
servations of detected gaps is advisable for a characterization of 
the planet. 

This work demonstrates that the gap carved by a moderately 
massive planet of one Jupiter mass at large orbital radii will be 
well within reach of ALMA. While the image quality naturally 
decreases with increasing distance, we find that the ability to ob- 
serve planet gaps is little affected by the source declination or 
inclination to the line-of-sight. Although few planets on wide 
orbits have been detected so far due to the biases of current tech- 
niques, their existence is supported by the more recent planet for- 
mation models as well as by observational constraints on planet 
populations. ALMA, probing a new part of the orbital param- 
eter space, has the potential to unveil more of them in nearby 
star-forming regions. 



Acknowledgements . We thank Frederic Gueth for enlightening discussions 
about simulating the effects of phase noise in ALMA reconstructed images and 
Remy Indebetouw, Scott Schnee and Juergen Ott for help on the use of the CASA 
ALMA simulator. This research was partially supported by the Programme 
National de Physique Stellaire, the Programme National de Planetologie and the 
Action Specifique ALMA of CNRS/INSU, France, and the Agence Nationale 
de la Recherche (ANR) of France through contracts ANR-07-BLAN-0221 
and ANR-2010-JCJC-0504-01. We acknowledge funding from the European 
Commission's 7 th Framework Program (contract PERG06-GA-2009-256513) 
and from the Millenium Science Initiative (ICM) of the Chilean ministry of 
Economy (Nucleus P10-022-F, "Disks with ALMA"). Simulations presented 
in this work were run at the Service Commun de Calcul Intensif (SCCI) de 
rObservatoi re de Grenoble, France. Images in Fig.^were made with SPLASH 
<Pricd2007l) . all other fi gures were made with the freeware Yorick. 



References 

Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 369, 216 
Alibert, Y., Mordasini, C, Benz, W., & Winisdoerffer, C. 2005, A&A, 434, 343 
Andrews, S. M., Wilner, D. J., Espaillat, C, et al. 201 1, ApJ, 732, 42 
Artymowicz, P. & Lubow, S. H. 1994, ApJ, 421, 651 

Barriere-Fouchet, L., Gonzalez, J.-E, Murray, J. R., Humble, R. J., & Maddison, 

S. T. 2005, A&A, 443, 185 
Bjorkman, J. E. & Wood, K. 2001, ApJ, 554, 615 
Boss, A. P. 2011, ApJ, 731, 74 

Burrows, C. J., Stapelfeldt, K. R., Watson, A. M., et al. 1996, ApJ, 473, 437 
Clampin, M., Krist, J. E„ Ardila, D. R., et al. 2003, AJ, 126, 385 
Cleeves, L. I., Bergin, E. A., Bethell, T. J., et al. 2011, ApJ, 743, L2 
Cossins, P., Lodato, G., & Testi, L. 2010, MNRAS, 407, 181 
Crida, A., Morbidelli, A., & Masset, F. 2006, Icarus, 181, 587 
Dodson-Robinson, S. E. & Salyk, C. 2011, ApJ, 738, 131 
Dutrey, A., Guilloteau, S., & Simon, M. 1994, A&A, 286, 149 
Fouchet, L., Gonzalez, J.-E, & Maddison, S. T. 2010, A&A, 518, A16 (Paper I) 
Fouchet, L., Maddison, S. T., Gonzalez, J.-F, & Murray, J. R. 2007, A&A, 474, 
1037 

Fukagawa, M., Hayashi, M., Tamura, M., et al. 2004, ApJ, 605, L53 
Gorti, U., Dullemond, C. P., & Hollenbach, D. 2009, ApJ, 705, 1237 
Grady, C. A., Polomski, E. E, Henning, T., et al. 2001, AJ, 122, 3396 
Grady, C. A., Woodgate, B., Bruhweiler, F. C, et al. 1999, ApJ, 523, L151 
Hashimoto, J., Tamura, M., Muto, T., et al. 2011, ApJ, 729, L17 
Hughes, A. M., Wilner, D. J., Calvet, N., et al. 2007, ApJ, 664, 536 
Lagrange, A.-M., Bonnefoy, M., Chauvin, G., et al. 2010, Science, 329, 57 
Lambrechts, M. & Johansen, A. 2012, A&A, 544, A32 
Lin, S.-Y., Ohashi, N., Lim, J., et al. 2006, ApJ, 645, 1297 
Lucy, L. B. 1999, A&A, 345, 211 

Marois, C, Macintosh, B., Barman, T., et al. 2008, Science, 322, 1348 
Marois, C, Zuckerman, B., Konopacky, Q. M., Macintosh, B., & Barman, T. 

2010, Nature, 468, 1080 
Mathis, J. S. & Whiffen, G. 1989, ApJ, 341, 808 
Mayor, M. & Queloz, D. 1995, Nature, 378, 355 

Mordasini, C, Alibert, Y., Benz, W., Klahr, H.. & Henning, T. 2012, A&A, 541, 
A97 

Muto, T., Grady, C. A., Hashimoto, J., et al. 2012, ApJ, 748, L22 



12 



J.-F. Gonzalez et al.: Planet gaps in the dust layer of 3D protoplanetary disks. II 

Otarola, A., Holdaway, M., Nyman, L.-A., Radford, S., & J., B. B. 2005, ALMA 
Memo, 512 

Owen, J. E., Ercolano, B., & Clarke, C. J. 2011, MNRAS, 412, 13 
Paardekooper, S.-J. & Mellema, G. 2004, A&A, 425, L9 
Pietu, V., Guilloteau, S., & Dutrey, A. 2005, A&A, 443, 945 
Pinilla, P., Birnstiel, T., Ricci, L., et al. 2012, A&A, 538, A114 
Pinte, C, Fouchet, L., Menard, F.. Gonzalez, J.-F, & Duchene, G. 2007, A&A, 
469, 963 

Pinte, C, Harries, T. J., Min, M., et al. 2009, A&A, 498, 967 

Pinte. C, Menard, F, Duchene, G.. & Bastien, P. 2006, A&A, 459, 797 

Price, D. J. 2007, PASA, 24, 159 

Quanz, S. P., Lafreniere, D., Meyer, M. R., Reggiani, M. M., & Buenzli, E. 2012, 
A&A, 541, A133 

Regaly, Z.. Juhasz, A., Sandor. Z., & Dullemond, C. P. 2012. MNRAS, 419, 
1701 

Schneider, J., Dedieu, C, Le Sidaner, P., Savalle, R., & Zolotukhin, I. 2011, 
A&A, 532, A79 

Varniere, P., Blackman, E. G., Frank, A., & Quillen, A. C. 2006, ApJ, 640, 1 1 10 

Wolf, S. & D' Angelo, G. 2005, ApJ, 619, 1 1 14 

Wolf, S., Gueth, E, Henning, T., & Kley, W. 2002, ApJ, 566, L97 



J.-F. Gonzalez et al.: Planet gaps in the dust layer of 3D protoplanetary disks. II, Online Material p 1 




5. 10.0 20. 50. 100.0 



Radius [AU] 

Fig. A.l. Map of the logarithm of the number of SPH particles 
per MCFOST grid cell in the meridian (R, z/R) plane for the 1 Mj 
planet. 
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Fig. A.l. Same as Fig. IA.ll in the horizontal (R, 0) plane. 



Appendix A: Dust density sampling in the MCFOST 
grid 

In order for MCFOST to produce accurate synthetic images from 
the output of SPH simulations, the dust density distribution must 
be sufficiently sampled with SPH particles in each cell of the ra- 
diative transfer grid. Fi gures lATI and IAT21 show 2D maps of the 
logarithm of the number of SPH particles per grid cell in the 
(R,z/R) and (R, <p) planes respectively. In most of the disk, the 
number of SPH particles contributing to the density in any grid 
cell varies from a few hundred to a few thousand (as a refer- 
ence, SPH codes typically use a few tens of neighbouring parti- 
cles to compute physical quantities). Only in the upper layers at 
large distances from the star does this number fall to values of a 
few, because there is very little material in these regions. Indeed, 
as can be seen in Fig. [TJ dust grains 100 ^m in size and larger 
settle efficiently to the midplane and smaller grains follow the 
distribution of the gas phase, which has a curved outer-rim and 
therefore a very low density in these regions. Both maps show 
smooth variations in the radial and vertical directions, whereas 
azimuthal variations are negligible (note the different colorscale 
in Figs. IA.ll and lA.2l ). and demonstrate the adequate sampling of 
the dust distribution for the radiative transfer calculations. 
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Fig. B.l. Radial temperature profile at the disk surface (solid red 
line) and in the disk midplane (black dashed line). 

Appendix B: Impact of the disk inner radius on the 
temperature structure 

The inner radius of the disk models is set to 4 AU in the SPH 
calculations, much larger than the characteristic inner radius of 
most T Tauri disks (usually located at the dust sublimation ra- 
dius as 0.05 AU). The resulting temperature structure is therefore 
incorrect in the central parts of the model, where the dust is di- 
rectly heated by the star. Our disk models remain optically thick 
in the radial direction however. As a consequence the temper- 
ature in the outer parts is computed correctly in these regions, 
as the heating of the disk midplane (where most of the millime- 
ter emission is coming from) is due to reprocessing of the stel- 
lar light absorbed and scattered by the disk surface. In these re- 
gions, the transfer of radiation is mostly vertical and the temper- 
ature structure does not depend on the details of the inner disk. 
Figure |B~T1 shows the radial profile of the surface and midplane 
temperatures. In the inner region, the dust is heated directly by 
the stellar radiation and both temperatures are equal. They start 
to decouple at r « 10AU, when the disk midplane no longer 
sees the star directly. The temperature structure and correspond- 
ing millimeter images can therefore be considered accurate from 
r > 15 AU, corresponding to 0. 1" in most of our simulations. 

Appendix C: Additional figures 

We present here figures similar to Fig. |4] for other wavelengths: 
A = 350 /im in Fig. IC.ll A = 1.3 mm in Fig. IC.2l and A = 2.7 mm 
in Fig.lOl 
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Fig. C.l. Same as Fig. H at A = 350 /mi. 
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